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1 .  INTRODUCTION 


Recent  (1990-1991)  military  operations  in  Southwest  Asia  resulted  in  the 
placement  of  materials  and  munitions  into  the  most  oppressive  combinations  of 
heat  and  humidity  in  the  world.  The  Persian  Gulf  is  a  shallow  body  of  water  less 
than  200  m  deep  that  attains  remarkably  high  temperatures  in  summer .  The  average 
water  temperature  in  July  and  August  is  31.1  °C  (88  “F)  .  One  of  the  highest  dew 
points  ever  recorded  was  33.9  °C  (93  °F)  at  Sharjah,  Saudi  Arabia,  which  is  near 
Dubai  on  the  southern  end  of  the  Persian  Gulf.  At  the  northern  end  of  the 
Persian  Gulf  exists  the  most  extreme  high  temperature,  high  humidity  environment 
on  record.  On  24  July  1953  at  1400  (LST) ,  the  temperature  was  48.3  °C  (119  "F) 
and  the  dew  point  was  29.4  °C  (85  °F)  ,  which  yields  a  relative  hiamidity  of  34 
percent . 

This  hostile  environment  could  adversely  affect  the  storage  life  and  operability 
of  military  equipment.  Among  the  requests  from  the  field  for  information  on 
environmental  conditions,  there  have  been  a  substantial  number  regarding  the 
interrelation  of  temperature,  pressure,  air  density,  and  relative  humidity.  This 
report  discusses  the  basic  equations  that  relate  these  variables  along  with 
simpler  approximations  chat  provide  the  accuracy  required  for  most  meteorological 
operations.  A  computer  code  that  computes  these  meteorological  parameters  is 
included  for  completeness.  This  computer  code  is  not  intended  to  be  used  for 
voluminous  data  input;  however,  modifications  could  be  made  so  chat  all  data 
could  be  input  in  one  routine  and  then  output  the  desired  results.  These 
equations  and  codes  were  supplied  to  AMSMC-QAN— P(D)  to  satisfy  an  immediate  need 
and  were  not  intended  for  field  use. 

2.  ESTIMATION  OF  HUMIDITY  PARAMETERS 

For  most  meteorological  applications,  the  moisture  content  of  the  atmosphere  is 
conveniently  measured  by  means  of  dry  (T)  and  wet  (T„)  bulb  temperatures . 
However,  the  computation  of  other  measures  of  humidity  such  as  the  dew  point 
(Tj) ,  vapor  pressure  (e) ,  and  relative  humidity  (RH)  can  be  somewhat  complex. 
Semiempirical  equations  (Goff  and  Gratch,  1945)  are  often  used,  but  for  many 
applications  simpler  and  less  accurate  equations  are  sufficient.  World 
Meteorological  Organization  Publication  8,  TP.  3  (WMO,  1969)  gives  the  universal 
psychrometric  formula  for  computation  of  humidity  parameters  as  the  Regnault 
equation  (equation  (1)). 


e  =  -  Ap  {T  ~  TJ 


(1) 


where  e  is  the  vapor  pressure,  p  is  air  pressure,  e„  is  the  saturation  vapor 
pressure  at  T„  (wet  bulb),  T  is  the  dry  bulb  temperature,  and  A  is  the 
psychrometric  coefficient.  The  saturation  vapor  pressure  is  the- source  of  the 
computational  problem.  The  difficulty  lies  in  the  departures  from  the  ideal  gas 
law  and  the  variation  of  latent  heat  with  temperature.  Several  papers  in  the 
literature  give  empirical  approximations  for  the  saturation  vapor  pressure 
(Murray,  1967,  and  Tabata,  1973,  for  example).  Saturation  vapor  pressure  can  be 
expressed  as  a  function  of  absolute  temperature  by  the  Claus is-Clapeyron  equation 
(equation  (2) ) . 
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1.  ^  =  -IL 

e,  dT  RT^ 


(2) 


where  e,  is  the  saturation  vapor  pressure,  T  is  the  absolute  temperature,  L  is 
the  latent  heat  of  vaporization  of  water,  R  is  the  gas  constant  for  dry  air,  and 
E  is  the  ratio  of  molecular  weight  of  water  vapor  to  that  of  dry  air. 
Integration  of  this  equation  is  difficult  because  L  is  not  constant,  but  varies 
with  temperature.  If  it  is  assumed  that  L  is  a  linear  function  of  T  then, 
according  to  Abbott  and  Tabony  (1985),  the  saturation  vapor  pressure  can  be 
expressed  as  in  equation  (3). 


e,  =  exp  (55. 17  -  6803  -  5 . 07  In  D 


(3) 


where  T  is  temperature.  From  Buck  (1981),  a  form  for  recovering  T  from  the  vapor 
pressure  over  water  is  given  in  equation  (4) . 


T  = 


CZ 

(B  -  Z) 


(4) 


where  Z  -  ln(e„/a)  and  a  -  6.1121,  B  -  17.368,  and  C  -  238.88  (over  ice;a  - 
6.1115,  B  -  22.452,  and  C  -  272.55).  Equations  (3)  and  (4)  appear  to  be  adequate 
for  most  meteorological  purposes.  Wexler  (1976)  optimized  these  equations,  and 
equation  (5)  is  a  most  accurate  formulation  for  e„,  with  uncertainty  being  only 
a  few  parts  per  million  (ppm) . 


=  0.01  exp  [-299 1.27  298"*  -  6017.0128e‘i  +  18.87643854 
-  0.0283  547210  +  0.17  838301  X  lO**  0^  -  0.84150417  X  10’’ 0^ 
♦  0.44412543  x  10'^*  0^  +  2.858487  ln0] 


where  9  is  absolute  temperature  in  kelvins.  Equation  (5)  also  provides  the  best 
available  estimate  for  vapor  pressure  over  supercooled  water  and  can  be  used  for 
water  values  below  0  *C.  Supercooled  conditions  are  often  encountered  with 
cooled  mirror  hygrometers  in  which  a  mirror  surface  can  maintain  a  layer  of  water 
at  -15  ®C.  The  equation  for  saturation  vapor  pressure  over  ice.  is  not  given  in 
the  text  but  is  included  in  the  computer  code  (see  appendix  for  the  "C"  program) 
for  completeness.  The  interested  reader  is  referred  to  the  paper  by  Wexler 
(1976)  for  the  exact  formulation. 

Sometimes  an  enhancement  factor  is  added  to  the  vapor  pressure  computation  to 
compensate  for  the  difference  in  saturation  vapor  pressure  over  "pure"  water  “to 
the  saturation  vapor  pressure  over  ambient  moist  air.  Hyland  (1975)  estimated 
this  uncertainty  at  1000  mbar  to  increase  from  0.01  percent  at  50  *C  to  0.06 
percent  at  -50  *C.  For  vapor  pressure  and  temperature  ranges  encountered  in 
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Southwest  Asia,  this  enhancement  factor  is  considered  superfluous.  The  interested 
reader  is  referred  to  a  paper  by  Buck  (1981)  if  the  factor  is  considered 
important.  For  most  meteorological  purposes  the  accuracy  required  of  e„  is  that 
equivalent  to  an  error  of  0.1  °C  in  temperature.  Equation  (5)  easily  satisfies 
this  requirement. 

3.  ESTIMATION  OF  THE  WET  BULB  FROM  THE  DRY  BULB  AND  THE  DEW  POINT 

The  precision  to  which  the  temperature  (T)  and  the  dew  point  (T^)  are  reported 
in  the  WMO  synoptic  code  messages  was  increased  in  1982  from  1  °C  to  0.1  °C. 
This  change  made  it  possible  to  compute  hiimidity  parameters  directly  from  the 
messages  rather  than  through  internal  reporting  of  the  wet  bulb  (T„)  .  The 
advantage  is  increased  automation;  but  the  recovery  of  T„  from  T  and  T^  is 
required,  and  this  procedure  is  not  straightforward.  Estimating  T„  from  T  and 
Tj  is  more  complex  than  obtaining  T^  from  T  and  T„.  Equation  (1)  is  e.xpressed 
in  terms  of  the  depression  of  the  wet  bulb  rather  than  the  dew  point,  and  so  can 
only  be  solved  for  T^  by  using  an  iterative  procedure.  A  simpler  method  is  to 
take  advantage  of  the  fact  that  T„  lies  between  T  and  T^.  At  low  temperatures 
T.„  is  not  much  less  than  T,  but  at  high  temperatures  T„  is  closer  to  T^.  Hence, 
the  ratio  of  the  wet  bulb  depression  to  the  dew-point  depression  always  lies 
between  zero  and  unity  and  increases  with  temperature.  Abbott  and  Tabony  (1985) 
give  a  convenient  linear  representation  of  this  property  in  equation  (6). 


T  - 


0.34  +  0.006  (r  +  Tj) 


(6) 


For  values  of  T  from  -10  ®C  to  +50  “C  and  dew-point  depressions  up  to  15  °C, 
errors  in  T„  from  using  equation  (6)  are  always  less  than  0.3  ®C.  For  increased 
accuracy  an  iterative  approximation  method  is  used  in  the  supplied  computer 
program. 

4.  ESTIMATION  OF  ABSOLUTE  HUMIDITY,  RELATIVE  HUMIDITY,  AND  AIR  DENSITY 

Calculation  of  these  values  begins  with  a  determination  of  e,,  which  requires  a 
knowledge  not  only  of  T  and  T„  but  also  of  e„  at  the  wet  bulb  temperature.  This 
can  be  obtained  by  replacing  T  by  T„  in  equation  (3)  to  yield  equation  (7). 


6.1070  exp 


17.38 
239.0  *  r„ 


(7) 


Since  e,  -  e  w^d)  •  Regnault  equation  may  be  written  as 


~  ~  Ap(T  -  T^) 


(8) 


The  ventilation  coefficient  A  is  usually  taken  to  be  0.00079  (0.000720  if  the  wet 
bulb  is  frozen).  Accurate  values  of  p  should  be  used.  Abbot  and  Tabony  (1985) 
show  errors  of  about  20  percent  in  e ,  2  percent  in  relative  humidity,  and  2  “C 
in  Tj  if  a  pressure  of  1000  mbar  is  used  when  in  reality  the  pressure  was  950 
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mbar.  The  correct  pressure  to  apply  is  the  site  pressure,  uncorrected  for 
altitude.  The  dew  point  can  be  calculated  from  the  transposed  Magnus  (1844) 
equation  for  evaporation  over  water  (equation  (9)). 


,  _  (239 .0  JO 

^  ~  (17 .38  -  JO 


(9) 


where  K  -  In  e  -  In  6.1070. 

The  saturation  vapor  pressure  at  the  dry  bulb  is  then  computed  so  that  the 
relative  humidity  (RH)  may  be  calculated  from  equation  (10). 


RH 


X  100 


(10) 


where  e^  is  the  saturation  vapor  pressures  at  T  and  Tj  . 

Absolute  humidity  (AH)  is  the  mass  of  water  vapor  per  unit  volume  of  atmosphere. 
The  equation  of  state  for  water  vapor  can  be  written  as 


( 

(RT,) 


(11) 


where  is  the  absolute  humidity  (grams  per  cubic  meter)  and  the  subscript 

refers  to  saturation  conditions  over  water,  e,  and  T,  are  the  vapor  pressure  and 
temperature  at  saturation  and  m^  is  the  mass  or  weight  of  the  water  vapor. 
Equation  (12)  is  a  convenient  computational  method  for  determining  absolute 
humidity . 

ah  =  1322^8314  Qjjp  [25.22  .y  (1.0  -  (273.15/rtf)  +  5.31  X  In  (273 . 15/rd)  J 


where  T  and  T^  are  air  temperature  (K)  and  dew-point  temperature  (K)  , 
respectively.  AH  is  expressed  in  grams  per  cubic  meter. 

Specific  humidity  (q)  is  defined  as  the  weight  of  water  vapor  contained  in  a  unit 
weight  of  air.  Units  of  q  are  usually  expressed  in  grams  per  gram  or  grams  per 
kilogram.  A  simple  expression  for  q  is  given  by  equation  (13). 


<7  = 


w 

1  ♦  w 


where:  w 


0.622e 
p  -  e 


(13) 
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where  e  is  the  vapor  pressure  ac  T  and  p  is  air  pressure  (mbar)  .  Units  are  grams 
per  gram. 

The  virtual  temperature  (T^)  is  defined  as  the  temperature  of  dry  air  having  the 
same  total  pressure  and  density  as  the  moist  air.  An  expression  for  T.^  is  given 
by  equation  (14).  p  and  e  have  the  same  meaning  as  in  equation  (13). 

r„=  r(l  +  0.61  w)  where-.  w=  °  ^  (1-4) 

p  -  e 


Once  T.,  is  known  the  air  density  can  be  calculated  from  equation  (15)  (List, 
1958)  . 


p  =  0.34838  ^ 


(15) 


where  p  is  air  density  in  kilograms  per  cubic  meter,  p  is  site  pressure  in 
millibars,  and  is  virtual  temperature  in  kelvins. 

5.  COMPUTATION  OF  PRESSURE  AND  DENSITY  ALTITUDES 

Pressure  altitude  (PA)  is  the  height  above  sea  level  in  the  standard  atmosphere 
at  which  the  altimeter  setting  (A)  exists.  The  PA  (in  meters)  can  be  computed 
by  using  equation  (16). 


PA  =  -^  [1  -  (A/P,)  He 


(16) 


where  T,  is  the  standard  sea  level  temperature  (288.16  K)  ,  a  is  the  standard 
lapse  rate  (0.0065  “’C/m)  ,  A  is  the  site  altimeter  reading  (inches  mercury  [Hg]), 
P,  is  the  standard  sea-level  pressure  (29.921  inches  Hg) ,  n  is  a  constant  with 
a  value  of  5.2561,  and  He  is  the  site  elevation  in  meters  above  sea  level. 
Restated,  pressure  altitude  is  the  distance  measured  from  the  29.921  inches  hg 
pressure  level,  which  is  the  standard  datum  plane. 

Density  altitude  (DA)  is  the  altitude  above  sea  level  characterized  by  a  known 
air  density.  Density  altitude  bears  the  same  relation  to  pressure  altitude  as 
true  altitude  does  to  indicated  altitude.  DA  may  be  computed  by  using  equation 

(17). 


DA  =  145,366  [1  -  17 . 326  {P,/T„) 


(17) 
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The  units  of  DA  are  in  feet;  and  A,  He,  P,,  and  n  are  defined  as  above  for 
pressure  altitude  and  P,,  is  defined  as 


=  [A*/"  -  /feCP,)*/"  (^)]'> 

^9 


6.  SUMMARY 

In  determining  the  various  measures  of  humidity,  the  greatest  difficulty  lies  in 
evaluating  the  saturation  vapor  pressure.  Thir  difficulty  is  related  to 
departures  from  the  ideal  gas  law  and  the  variation  of  latent  heat  with 
temperature.  These  departures  make  the  Claus ius-Clapeyr on  equation  difficult  to 
integrate  analytically.  Semiempirical  equations  by  Goff  and  Gratch  (1946)  are 
available  for  extremely  accurate  computations  of  saturation  vapor  pressure,  but 
for  most  meteorological  purposes  this  accuracy  is  superfluous.  Equations  given 
in  this  report  are  simpler  approximations  providing  the  accuracy  required  for 
most  meteorological  applications.  A  simple  iteration  routine  has  been  supplied 
for  determining  the  wet  bulb  temperature  from  the  temperature  and  dew  point. 
From  the  equations  and  the  supplied  computer  program,  most  vapor  pressure, 
humidity,  air  density,  and  altimetry  computations  can  be  accomplished.  P-sults 
compare  favorably  with  values  in  meteorological  handbooks. 
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APPENDIX 


"C"  PROGRAM 


File:  S7P.C 

Pro^raa:  Saturation  Vapor  Pressure. 

Author:  Kevin  X.  Lems.  1991. 

Purpose:  SVF  computes  the  saturation  vapor  pressure  of  a  set  of  input 

data.  See  the  provided  documentation  for  a  full  explanation. 

SVP  was  written  for  the  Weather  Data  Indegration  Branch  of  the  Battle 
Weather  Data  Division  of  the  Atmospheric  Sciences  Laboratory. 

S7F.C  was  compiled  using  Borland's  Turbo  compiler.  It  was  not 
compiled  using  any  options,  so  it  should  compile  with  any  Borland 
compiler  and  run  properly. 

STP.C  should  easily  port  to  any  system  with  an  AXSI  C  compatible 
compiler. 


tinolude  <3tdio.h> 
ilnclude  (math.h) 

void  Svp(  double  *7,  double  >Tv,  double  'P  ); 
double  Sw(  double  T  ) ; 
void  Alr_Density(  double  T.  double  P  ); 
void  Tempo ; 

void  Wet.Tf  double  T,  double  Td,  double  P  ); 
double  Dem.Polntf  double  T,  double  P  ); 
void  Density^Altltudef  double  iT  ); 
void  Pressure.AltitudeO ; 

Int  Qet.Dataf  double  >data,  cbar  param,  cbar  tt  ); 
double  Convert  .Temp  (  double  temp,  char  type,  char  need  ); 
double  Convert_Pressure (  double  p,  char  type,  char  need  ); 
double  Convcrt,Beight(  double  h,  char  type,  char  need  ); 
char  QetcharO; 

/>  Define  the  varaibles.  >/ 
double  X,  ewp,  Td,  Tv,  B,  rbo; 
char  numatrCSO],  ch; 


void  aainO 
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double  T.  Tv.  P: 
pute( 
puts! 
puts! 

puts!  ■''nInstructions:\n' 


»»»  Saturation  ?apor  Pressure  »»» 


1; 


puts!  'When  prompted,  lollow  your  input  by  a  letter  corresponding  to  the  units  of  ); 
puts!  ‘the  input  values.  For  example,  20  degrees  Celsius  is  entered  as  20c  (capital' 
puts(  ‘or  lower  case  'O.  293  Kelvin  as  293X.  68  degrees  Fahrenheit  as  dSF.  29.92'  ); 
puts(  ‘inches  of  mercury  as  29.921,  1013  millibars  as  1013a,  1000  feet  as  lOOOF,'  ); 
puts!  'and  300  meters  as  300li.  Output  units  are  given  nth  the  computed  value.'  ); 


do  ( 


Svp(  iT.  Hv,  iP  ): 

Air_Density(  Tv,  P  ); 

do  { 

print! (  '\nWould  you  like  to  compute  the  Density  Altitude  (y/n)'>  '  ); 

ch=Qetchar() ; 

if(  oh!=’y’  a  chi-’Y'  A&  chi=’n’  t&  chl=’H’  )  ( 
putchl  '\a'  ); 
chsO; 

) 

)  while!  'ch  ) : 

if (  ch"'y’  1 !  ch*=’T'  ) 

Densily_iltitude{  H  ); 

do  ( 

prlntf!  ‘\nlfould  you  like  to  compute  the  Pressure  Altitude  (y/n)'  '  ); 

ch^OetcharO ; 

lf(  chl=’y'  U  chl-’Y'  U  chl^’n'  Al  ehl-’i’  )  { 
putch(  ’\a'  ); 
ch=0; 


)  while!  Ich  ) ; 
if!  ch«’y’  1!  ch=='Y'  ) 
Pressure_Altltude!) ; 


...  We  are  done  with  all  the  computations.  Find  out  if  the  user  has  some 
other  stuff  to  do. 

do  ( 

printf!  '\n\nWould  you  like  to  co^iute  additional  values  (y/n)?  '  ); 
eh=Oetchar ! ) ; 

if!  chl=’y’  U  chl=’T’  U  chl=’n’  U  chl=’I’  )  ! 
putch!  '\a'  ); 
ch=0: 

) 

}  while!  Ich  ) ; 

)  while!  chss’y’  1 :  ch*s’Y'  ) ; 


} 
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void  Svp(  double  »T,  double  »Tv,  double  »P  ) 


( 

double  Afl,  Hfl,  q,  »,  e; 

/«  Get  the  dry  bulb  teaperature.  ♦/ 
puta(  'Xn'  ); 
do  ( 

printK  "Input  the  dry  bulb  teaperature  T  followed  by  K,  C,  or  F:  '  ); 

)  while!  Qet_Data(  T.  ich  )”-l  ); 

»T-Convert_Teap(  »T,  ch,  ’K’  ); 

/*  Get  the  pressure .  »/ 
do  ( 

printK  "Input  the  pressure  P  followed  by  I  or  M  (default=1013.2Sab) ;  "  ); 

)  while!  Get_Data!  P,  Ich  )==-!  ); 

»P=Convert_Pres3ure!  »P,  ch.  'll’  ); 

do  { 

prlntf!  "\nDo  you  have  the  dew  point  teaperature  !y/n)?  "  ); 
ch=Getcbarl) ; 

IK  chl  =  ’y’  ii  chi-’Y’  ii  chl-'n'  ii  cbl-’V  )  ( 
putch!  ’\a'  ); 
eh=0: 

) 

}  while!  'oh  ) : 
putchar!  10  ); 
if!  ohs-’y'  II  ch=='T’  )  ( 
do  ( 

printK  "Input  the  dew  point  teaperature  Td  followed  by  X,  C,  or  F:  "  ); 

}  while!  Get.Datal  Md,  'T',  ich  )=s-l  ); 

Td^ConverWeap!  Td,  ch,  ’X’  ); 
e*Ew(  Td  ) : 

Wet  Tl  «T,  Td,  IP  ); 

) 

else  ( 

do  ( 

printK  "Input  the  wet  bulb  teaperature  Tw  followed  by  X,  C,  or  F:  "  ); 

)  while!  Oet.Data!  iTw,  'T',  ich  )-=-!  ); 

Tw=Convert.Teap|  Tw,  ch,  'X'  ); 

Td=Dew_Point(  iT,  iP  ); 

prlntf!  "\nlhe  Dew  Point  Teaperature  is  \t\tX10.4f  [Xelvln]\n',  Td  ); 

\t\tX10.4f  [Celclus]\n' ,  Convert  Tewp!  Td,  'X',  'C  ) 

'  \t\tX10.4f  [Fairenhelt]\n" ,  Convert  Teap!  Td  ’X’  ’F’ 

e^8w(  Td  ) :  -  -T  . 

) 

/>  Coapute  'cap'  •/ 
ewp=gw(  iT  ); 

prlntf!  "\nlhe  Saturation  Yapor  Pressure  e!T)  is  \tX10.4f  [illllbarslXn",  ewp  ); 
prlntf!  "The  Saturation  Vapor  Pressure  e!Td)  is  \tX10.4f  [alllibars]\n" ,  e  ); 
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/>  Coipute  the  relative  huaidity.  •/ 

HH=e/eiip»100: 

printfl  'The  Relative  Huaidity  is  \t\t210.411X\n‘,  SH  ); 

/<  Coapute  the  absolute  huaidity.  >/ 

iH=(  1322.8314/(»T)  )»exp(  25.22M  1.0  -  (  273.15/Td  )  )  +  5.31»log(  273.15/Td  )  ); 
print!!  ’NnThe  Absolute  Huaidity  is  \t\lI10.4l  l^raas  per  cubic  aeter]\n’,  AH  ); 

/*  Coapute  the  virtual  teiperature.  >/ 
f=0.822»eap/(  »P  -  eap  ) : 
q=*/(  1.0  ♦  «  1 : 

print!!  'The  Speelllo  Huaidity  Is  \t\t210.4!  (Iraas  per  graa]\n‘,  q  ); 

«Tv=»T»!  1.0  ♦  0.61»a  ) ; 

print!!  'The  71rtual  Teaperature  is  \t\tll0.4!  (Kelvin]\n‘.  »Tv  ); 


double  £«!  double  T  )  /•  T  aust  be  in  K,  and  P  aust  be  in  aB.  >/ 

( 

double  satvp; 

/I 

The  second  equation  Is  used  for  saturation  vapor  pressure  over  Ice.  I!  you 
wish  to  use  this  equation,  remove  the  coaaent  aarkers  on  the  indicated  lines. 

1/ 

/i  Remove  this  line, 
i!!  T)=273,1S  ) 

Beaove  this  line.  •/ 

satvp--0.0Uexp!  -2991.2729/!  TiT  )  -  6017.0128/T  ♦  18.87643854  - 
0.028354721»T  ♦  0. 17838301e-4»T»T  -  0.84150417e-9iTiT»T  ♦ 
0.44412543e-12iTiTtTiT  ♦  2.858487iloj(  T  )  ): 

/•  Remove  this  line, 
else 

satvp=0.01»exp(  -5865.3696/7  ♦  22.241033  ♦  0.013749042iT  - 

0.34031775e-4iT»T  ♦  0.269e7687e-7iTsT»T  ♦  0.8918651ilo8 I  T  )  ); 

Remove  this  line.  >/ 

return  satvp; 


/I  Coapute  the  air  density.  >/ 
void  Alp_Denslty!  double  Tv,  double  P  ) 

( 

double  rbo; 
pho=0.34838»P/Tv: 

print!!  ‘\n7he  Air  Density  is  \t\t\tZ10.4!  (kllojraas  per  cubic  metcrlVn',  rho  ); 


) 
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it  Conpute  the  dew  point.  »/ 
double  Dew-Point!  double  T,  double  P  ) 

( 

double  A=0. 000799,  es,  e.  X; 

T  -=  273.18;  Tw  -=  273.16;  /»  convert  to  deg  C  »/ 

es=a.l07»exp(  17.38»7w/(  239.0  ♦  Tw  )  1; 
e  :  es  -  A»P* (  T  -  Tw  ) ; 

X  =  log(  e  )  -  log(  8.1070  ); 

/»  printfC\ne,  es,  X.  P,  T.  Tw  =  17.31,  X7.3f,  ;7.3f.  Z8.3f.  17. 3f.  X7.3f\n', 
e,  es.  X,  P,  T.  Tw) ;  »/ 
return  (  239. 0*X  )/(  17.38  -  X  )  ♦  273.16; 


} 


void  iltet_T(  double  T,  double  Td,  double  ?  ) 


Int  tlaee=0; 

double  twO,  F,  Fp,  dlf; 

double  A,  a,  b.  c; 

T=Convert_Teiip(  T,  'X',  ’C  ); 

Td=Convert.Teap(  Td,  ’X’,  'C  ); 
tw0s0.3»(  T  ♦  Td  ); 
do  ( 

IK  tw0>=0.0  )  { 

A-0. 000799; 

4=17,38; 

b»239.0: 

) 

else  ( 

A=0. 00072 : 

4=22. 44; 
b=272.4; 

) 

c=8.107«exp(  a»Td/(  b  ♦  Td  )  )  ♦  A»P»T; 

F  =  8.107iexp(  a«tw0/(  b  ♦  twO  )  )  ♦  A«P»twO  -  c; 

Fp  =  8.107ia»b/(  (  b  ♦  twO  )»(  b  ♦  twO  )  )iexp(  a»twO/(  b  ♦  tig)  )  )  ♦  AiP; 

dif=F/Fp: 

do  ( 

dK=0.5»(iif; 

}  wbile(fabs(dlf/tw0)>=0.5); 
two  =  two  -  dll; 

IK  tiMsee)20  )  ( 

puts(  ’\niarnlng:  possible  error  In  co^iutlng  wet  bulb  ten|>erature"20  Iterations  done.*  ); 
breaJi; 

) 

)  while!  labs!  dll  )>=0.001  );  /•  11  It's  accurate  enough;  stop.  «/ 

prlntK  ‘\nThe  wet  bulb  te^ierature  Is  \t\tX10.41  (CelsiuslNn' ,  twO  ); 

prlntK  ■  \t\tX10.41  [Celslusl\n’,  Convert_Teip(  twO,  ’C’,  ’X’  1  ); 

printK  ■  \t\tX10.41  (Fahrenhelt]\n' ,  Conwert_Teip(  tirt,  'C,  'F'  )  ); 

7w=tw0: 


17 


) 


void  Density_iltitude(  double  •?  ) 


double  rvOOS,  a=.0065,  n-S.25ai,  Ps=29.921.  T3-288.18; 
double  Alt,  Be,  rbo_h,  Pm,  Tv; 
cbar  numstrlSOl; 

pplntK  'Nnlnput  T  lollomed  by  K,  C,  op  F  (default  is  oplglnal  lnput=I7.2{K) :  ’,  »T  ) ; 

IK  Qetjatal  T,  'T',  ieh  )l=-l  ) 

*T=Convept_TeBp(  »T,  cb,  'K'  ); 

Tv  =  -0.288  »*(  9. 0/5.0  )»(»T)«(  1.0  ♦  1.60779ip  )/(  1.0  ♦  p  ); 
do  ( 

ppintK  ‘Input  Site  Pressure  Oncorrected  follomed  by  I  op  M:  ‘  )  ; 

)  rtilel  Getjatal  lAlt,  'P',  ich  )==-!  ); 

Alt=Convept_Ppessure(  Alt,  ch,  T  ); 
do  { 

ppintK  'Input  site  elevation  (ollomed  by  H  or  F:  ‘  ); 

1  mhilel  Get_Data(  Me,  'H',  Ach  )==-!  ); 

He=Convept_Height(  He,  ch.  'M'  ); 

Pm=pow(  pow(  Alt,  1.0/n  )  -  8e*po»(  Ps,  1.0/n  )•(  a/Te  ),  n  ); 
pho>  145388. 0»(  1.0  -  pow(  17.328»(  Pm/Tv  ).  0.235  )  ); 
pplntK  ’The  Denelty  Altitude  ii  \t\tX10.4f  [feet)\n‘,  pho,h  ); 

pplntK  ‘  \t\tX10.4f  (ieteps]\n‘.  Convept_Heljht(  phoh,  ’F’,  ’M’  )  ) 


void  Ppeuure.AltitudeO 


double  PA,  Ts,  a.  A,  Ps,  n.  He; 

T8=288.18; 
a=0.0085: 
puts!  '\n'  ); 
do  ( 

pplntK  ‘Input  the  Altlictep  Value  (olloved  by  I  op  M:  '  ); 

)  mhllel  OetJaUl  U,  ’P',  &ch  )==-!  ); 

A=Convept_Ppessupe(  A,  ch,  T  ); 
do  ( 

pplntK  'Input  the  site  elevation  folloacd  by  M  op  F;  '  ); 

)  ifiile(  Get.Datal  Me,  'H',  Ach  )»-l  ); 

Be^Convept. Height!  He,  ch,  'M'  ); 

Ps*29.931; 

n*5.25fll; 

PAK  Ts/a  )*(  1.0  -  pow(  A/Ps,  1.0/n  )  )  »  He; 

printfl  ‘The  Pressure  Altitude  is  \t\tZ10.4f  [aeterslNn',  PA  ); 

printK  '  \t\tZi0.tf  (feet]\n‘.  Convert.Height(  PA,  'M',  7'  )  ); 


) 
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int  3et,Data(  double  (data,  char  paras,  char  'type  ) 


{ 

int  len; 
double  tip; 
cbar  t[50]; 
jetei  t  1 ; 

ill  strempt  t,  "  )==0  )  ( 

Iff  parai==’P’  )  ( 
idata=1013.2S; 
itype=’ll’ : 
return  0; 

) 

elae 

return  -1: 

) 

len=gtrlen(  I  ); 

»type=t[len-ll: 
switch!  parai  )  ( 
case  ’7’: 

if!  »typel  =  'C’  SA  »typel=:'e’  H  (typei  =  ’X’  &&  »type'  =  ’k'  itype'  =  ’F’  &&  (type':’!’  ) 
return  -1; 
break; 
case  'P': 

if!  itypeis'li’  ii  (type's’!’  li  (type's’I’  ii  (type's’!’  ) 
return  -1; 
break; 
case  ’H*; 

if!  (typeis'M'  U  (typels'i'  ik  (typels’F’  U  (typels’f’  ) 
return  -1; 
break; 

) 

t(len-l]s’\0’: 

(datasatof It); 

tapsConvert.Teip!  (data,  (type.  'C  ); 
if!  paraiss’T’  )  { 

if!  t^)<-ioo.o  ::  tipiioo.o  )  ( 

puts!  'Te^ierature  ia  out  at  ranje.  Enter  a  new  teiperature,  -100'C(T! lOO'C  ); 
return  -1; 


) 

return  0; 


) 


double  Convert  .Teip  !  double  teip,  cbar  type,  cbar  need  ) 
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( 

switcbl  need  )  ( 
cue  'K': 

switch (  type  )  ( 
case  ’k’: 
case  'K’: 

return  te^); 
case  ' f ' ; 
case  'F* : 

return  5.0/9.0»(  leap  *  32.0  )  ♦  273.15; 
case  'c': 
case  'C: 

return  tewp  *  273.15; 

) 

case  'C ; 

switch!  type  )  ( 
case  'k': 
case  'X': 

return  tewp  -  273.15; 
case 

case  'F': 

return  5.0/9.0»(  tewp  -  32.0  ); 
case  'c' ; 
case  'C; 

return  tewp; 


case  7'; 

switch!  type  )  ( 
cue  'k'; 

cue 

return  9.0/5.0*(  tewp  *  273.15  )  *  32.0; 
case  ’c’: 
case  'C : 

return  9.0/5.0*tewp  *  32.0; 
case  ' f ' : 
case  7*; 
return  tewp; 

) 

) 

return  -1; 


) 
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double  Convert, PTessurei  double  p.  char  type,  char  need  ) 


( 

switch!  need  )  ( 
caae  ’ I ' : 

switch!  type  I  { 
caae  ’ 1’ : 
case  T: 

return  p; 
case  ’b' : 
case  'M’: 

return  p/33. 86; 

) 

case  'M': 

switch!  type  )  ! 
case  '  i ' : 
caae  T : 

return  33.86»p; 
caae 

case  'M' ; 
return  p; 

) 

) 

return  -1; 


1 


double  Convert, Height!  double  h,  char  type,  char  need  ) 

! 

switch!  need  )  ( 

caae  'M'; 

switch!  type  )  ! 

cue 

case  7': 

return  b*0.3048; 

cue 

cue  'M'; 
return  b; 

) 

cue  7': 

switch!  type  )  ( 

cue 

cue  7’: 

return  h; 
cue  'i' ; 

cue  'MV 

return  h*3.28083080S; 

) 

) 

return  -1; 
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t 


char  QetcharO 

( 

char  o: 
c^getcharO : 
ffluahi  stdln  ) 
return  o: 

} 
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